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i ABSTRACT 

bX) 1 Basing on the fundamental symmetry that the space-time inversion is 

^ | equavalent to particle-antiparticle transformation, a relativistic modifi- 

<^ ■ cation on the stationary Schrodinger equation for many-particle system 

is made. The eigenvalue in the center of mass system is no longer equal 
^ ■ to the negative of binding energy simply. The possible applications in 

various fields (e.g. the model of quarkonium) are discussed. 
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As is well known in quantum mechanics, for a system composed of two particles with 
mass mi and m 2 and interaction potential V(\f[ — | ) , (f[ and r~2 being the coordinates 
in laboratory (L) system), after introducing the coordinate of center of mass (CM) R = 
jf(mif[ + m 2 r2), (M = mi + m 2 ), and relative coordinate r = — f[, one easily obtain 
the stationary Schrodinger equation 

[-^Vl+V{r)}4>{r)=e4>{r) (1) 

with reduced mass \i = m^/M and eigenvalue 

e = E- — -Mc 2 (2) 

— * 

where E is the total energy of system in L system and P = \P\ is the momentum of CM. 
When P = 0, e = — B < 0, B is called as the binding energy of system. 

Usually, it is said that the above "nonrelativistic approximation" is good for the case 
of P 2 /{2M) being small. However, this is doubtful. One can always set the laboratory 
into motion to render P arbitrary large. Then the accuracy of calculation in CM system 
would rely on the motion of external world ! Does it respect to the "principle of relativity" 
in special relativity (SR) ? The reason why one can not isolate his calculation from the 
external world is the following. The summation of e with P 2 / (2M) is made along a 
straight line in Eq.(2) whereas the correct relation in SR is that of a right triangle: 

E 2 = P 2 c 2 + (Mc 2 - B) 2 (3) 

Hence, the statement that "the total energy E is equal to the sum of kinetic energy 
(P 2 /2M) of CM and the internal energy (Mc 2 — B)" is not rigorous. One should say that 
"E is equal to the square root of the sum of square of kinetic energy (Pc) of CM and 
that of internal energy ". Now the problem is "how to modify the Eq.(l) for meeting the 
requirement of Eq.(3) ?" 

It is said in ancient China that "to gain new insights through restudying old material" . 
Let us restudy an alternative derivation of Klein-Gordon (KG) equation. Consider a 
spinless particle with rest mass mo being in free motion. Then its wave function 9(x,t) is 
described by a "nonrelativistic quantum equation" : 

d k 2 
ih—6(x, t) = m cl9(x, t) - - — V 2 6(x , t) (4) 
at IrriQ 

Here, we add a term of rest energy moc 2 with c\ being merely a (unfixed yet) constant 
with dimension of velocity. The next crucial step is assuming that inside a particle state 
9, there is always a hiding antiparticle state £). 9 and x are coupled together via 
motion. Instead of Eq.(4), we should have a simultaneous equation as follows: 

l h§ l 9 = moc?* - £^V 2 9 - £-V 2 x 

^§- t X = -m c 2 lX +£;V 2 X+^ 2 e [ ) 



The guiding rule for establishing the Eq. (5) is a basic symmetry. It is invariant under the 
space-time inversion x — > —x, t — > —t and transformation: 

x (x,t) = 9(-x,-t) (6) 

For solving Eq.(5) in general, we use the ansatz 



-< 



= U + i-K- 



(7) 



and obtain the K-G equation 



Its plane wave solution 
leads directly to 



^- v2 + ! f i) ^ i) ^ (8) 



(f>(x, t) = exp[i(p ■ x — Et)/h] (9) 



E 2 = ?<* + m 2 ci (10) 

For clarifying the meaning of Ci, we look at the velocity (v) of particle which is equal to 
the group velocity (v g ) of de'Broglie wave: 

where the quantum relations E — ftuj and p = hk have been used. 
The inertial mass is defined as 

p ,,dE. 1 d n 

"* = H /( * ) = 23^ (12) 

Combining Eqs.(10)-(12), we arrive at 

E = mcj (13) 



v 2 



m = m /\l j 



c 



as expected. But here c has the meaning as the limiting speed of particle, its value is 
obtained from the measurement on the n meson beam and is coinciding with the speed 
of light, c: 

Cl = c = 3 x 10 10 cm/sec (15) 

In the above derivation of K-G equation (a reversed version of that in Ref.[l]), we start 
from E\ v= o = m c^ and get eventually the mass energy relation E = mc\. The proof bears 
some resemblance to the "inductive method" in mathematics. However, the important 
thing is injecting into the proof a "relativistic principle", i.e., the basic symmetry (6), 
which plays the role of "hormone" for activation of mass from tuq into m. 



The symmetry (6) is discussed generally as a statement that "the space-time inversion 
is equavalent to particle-antiparticle transformation" in Refs. [2-3]. In our opinion, it 
is a natrual postulate after we learn carefully from the development of physics since the 
discovery of parity violation [4,5] and the observation of Schwinger et al. [6,7]. 

We are now in a position to generalize the above derivation to two particle CclSC clS 
in Eq.(l). Denoting 9 = 9(f[,f 2 ,t) and x — xifii^2,t) the particle and correspond- 
ing "antiparticle" state again, then instead of Eq.(5), we can write down the following 
simultaneous equation: 

ih% = Mc*9 - (^V 2 , + £p%){6 + X ) + V(\fl - r 2 \)(9 + x ) 
in^^-Mc\ + (£-Vj. 1 + £-Wj,)(9 + X )-V(\r 1 -r 2 \)(0 + x) 1 J 

which still respect to the symmetry: 

X(r~i,f 2 ,t) = 9(-f{, -f 2 , -t) (17) 

As Eq.(7), we set 

" = *+ i s^*. x-*-^* < 18 > 

with $(fl,f2,t) — > $(R,r,t) obeying the equation: 

1 

— V^> + - 2 

n n 

Factorizing the solution as 

$(R, f, t) = e ip - R ' h e- iEt / n ^{r) (20) 
and substituting it into Eq.(19), we arrive at 



$ - c 2 V|$ - c 2 — V 2 $ + (M 2 c 4 + 2VMc 2 )§ = (19) 



l-iVl+V(r)]m = em 
e= 2 jh(E*-M^-P^) 



(21) 



Note that the eigenvalue e ^ E — Mc 2 even when P = 0. Comparing Eq.(8) with 
Eq.(3), we find that the accurate relation between the binding energy B and e reads 

5 = Mc 2 [l-(l + ^)^] (22) 

Only when << 1, can one recover the nonrelativistic approximation: 

B ~ -e (23) 

In general case, one should use Eq.(22) or its equavalent (for P = 0): 

E\ P=0 = [2Mc 2 e + M 2 c 4 ] 1/2 (24) 

It is easy to generalize the above consideration to many particle (n > 3) case. Denote 
the coordinates of ith particle with rest mass m 8 are r\ and r*j in L system and CM 



1 

i v^ivi itscius i\, — 2- 

while r*j = r'i — R obeys the constraint: 



system respectively. The coordinate of CM reads R = J2?=i ~tr m i r 'ii (M = Y!i=i m i)i 



E m ^ = ° ( 25 ) 
i=i 

Direct calculation leads to 

2=1 £ 2=1 1 2=1 J — 1 

The third term in RHS can be discarded because in CM system the total momentum 
equals to zero in stationary state. Denoting 9 = 9{r\, r' 2 , • • • , r' n ), x = x( r 'ii r 'i-, • • • , r' n ) 
and introducing again the wave function 

$(R, f u t) == e*™^-™'*^, ? 2 , ■ ■ ■ , r n ) (27) 

we find 

J [-¥ sr=i + = ei, 

Similar to Eq.(21). However, only (n — 1) coordinates in r*j are independent due to the 
constraint (25). 

Some remarks are in order: 

1. There is something unreasonable in previous stationary Schrodinger equation (1). The 
eigenvalue e may fall downward without lower bound. This situation does occur for a 
singular potential like V(r) ~ ^ (see Refs.[8-9]). Now Eq.(21) has no this kind of worry. 
There is a minimum value for e: e min = — |Mc 2 , or E^ n = 0. Therefore, the solution 
of Eq.(21) may be viewed as a variational problem and a lower bound exists for any 
variational procedure. 

2. Let us compare Eq.(21) with Dirac equation. For an electron with mass m e = m 
moving in the Coulomb field V(r) = —Ze 2 /(Airr) of neucleus with mass itln — > oo, the 
total energy of electron reads [10]: 

Ed = mc 2 {\ + Z2a2 = Y 1 ' 2 (29) 

[/(j + l/2) 2 -2V + n'f 

with a = ~ j|=, j = |, |, • • •, n' = 0, 1, 2, • • •. In terms of the principal quantum 



number n — n + (j + |), one has 



,-, l(Za) 2 l(Za) 4 1 3 



E D = mc 2 \l - r - "H 1 ( 30 ) 

[ 2 n 2 2 n 3 v j + l/2 4n ; J v ; 

For comparision, the rest energy of neucleus tonc 2 must be substracted from the E derived 
from Eq.(24). Denote 

E S = E- m N c 2 = c 2 {m + m N )[l + 26 . J 1/2 - m N c 2 (31) 

(to + ton)c 



where e is well known as 



#s = mc\l - * mN A * - * rnm% Za 4 _ _ 

1 2(m + mjv) V n ; 8(m + mjv) 3V n ; J V ; 

The difference between (30) and (32) is stemming from two reasons. In Dirac equation 
the spin of electron is taken into account whereas in Eq.(21) the finiteness of neucleas 
mass is important. However, both equations are relativistic because the basic symmetry 
(6) is respected (see Refs.[2-3]). 

3. Being an improvement to Eq.(l), Eq.(21) or (28) brings some modification on many 
problems in stationary states. Roughly speaking, for binding state problem, the modifi- 
cation is very small in atomic physics (< 10 -5 ), it accounts for 10 -3 in nuclear physics 
and may reach 10~ 2 in particle physics. For the high energy scattering problem, this 
modification might be also important. 

4. For concreteness, let us have a quick look at the potential model of heavy quarkonium, 
QQ. Assume that the potential between Q and Q is of the linear type, V(r) = ar, with 
constant a being independent of quark mass m. For S states, the stationary equation is 
solved analytically with eigenvalue 

en = A„(|V 3 , (n = l,2,...) (33) 

A n being the zeros of Airy function [11]. So the energy (mass) of QQ reads from Eq.(24): 

K = 4Ml + y(|i) 1/3 ] 1/2 04) 
with fj, — y — ^. On the other hand, the previous equation (1) with (2) yields 

K = V + A„(|> 3 (35) 

Table 1 is a comparision between the measured energy E^ p of S states in Upsilon bb [12] 
and E n or E n . In either case, the parameters a and \i are adjusted so that E n or E n is 
coinciding with E^ p for n — 1 and 2. the mass of constituent quark b is fitted as 

m b = 2fx = 4.326GeV (36) 

a = 0.4530GeV 2 

from E n versus 

m ' b = 2/i' = 4.354GeV (37) 
a = 0.3804GeV 2 



from E' n . The general trend of E n for higher n seems better than that of E' n as shown in 
the table. 

5. Similar fitting procedure used for Charmonium J/ip = cc by Eq.(34) leads to 

ra c = 1.031GeV (38) 
a = 0.4183GeV 2 

whereas Eq.(35) yields 

m' c = 1.155GeV (39) 
a = 0.2099GeV 2 

6. If neglecting the dependence of constant a in V(r) = ar on the quark mass, we may 
discuss the dependence of quarkonium mass on the quark mass for the level with same 
quantum numbers. The Feynman-Hellmann theorem for stationary Schiodinger equation 
reads [13] 

! = -I(e-<V))<0 (40) 

Now it is replaced by 

dE 4 / /,m 8// E 

The latter two terms in RHS are positive. Combining Eq.(41) further with the virial 
theorem: 

i ,71 ' 

(42) 

(43) 
(44) 



we get for V(r) = ar: 
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from E' = e + 4/i'. Eq.(44) can easily be integrated as 

E'(fi') =4//' + CV~ 1/3 (45) 
whereas Eq.(43) can also be linearized by E = ^fy and integrated as 

E(fj) = [16/i 2 + C/i 2/3 ] 1/2 (46) 

As an interesting test, we use the experimental data of ground state energy for Charmonim 
E 1 = 3.097 GeV to fix the value of C in Eq.(46) or C in Eq.(45) (using the value of fx or 
fi in the Eqs.(38) or (39) at the same time). Then the ground state of bb can be estimated 

as 

f?i(66) = 9.420GeV (47) 

or 

^(66) = 9.214GeV (48) 



versus 

E™ p (bb) = 9.460GeV (49) 

Similarily for the 2S state (3.686 GeV) of Charmonium with other value of C or C', we 

get 

E 2 (bb) = 9.9562GeV (50) 
E' 2 (bb) = 9.5922GeV (51) 

versus 

E? p {bb) = 10.023GeV (52) 
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Table 1: The S states of Upsilon bb. 



n 


1 


2 


3 


4 


5 


6 


£ n xp (GeV) 


9.46037 


10.023 


10.355 


10.580 


10.865 


11.019 


E n (GeV) 


9.46037 


10.023 


10.461 


10.834 


11.163 


11.462 


E' n (GeV) 


9.46037 


10.023 


10.483 


10.890 


11.262 


11.609 


A n 


2.338 


4.088 


5.521 


6.787 


7.944 


9.023 
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